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Abstract. Magnetic clouds (MCs) are the interplanetary coun¬ 
terparts of coronal mass ejections (CMEs), and usually modeled 
by a flux rope. By assuming the quasi-steady evolution and self¬ 
similar expansion, we introduce three types of global motion into 
a cylindrical force-free flux rope model, and developed a new 
velocity-modifled model for MCs. The three types of the global 
motion are the linear propagating motion away from the Sun, 
the expanding and the poloidal motion with respect to the axis 
of the MC. The model is applied to 72 MCs observed by Wind 
spacecraft to investigate the properties of the plasma motion of 
MCs. First, we And that some MCs had a signiflcant propagation 
velocity perpendicular to the radial direction, suggesting the di¬ 
rect evidence of the CME’s deflected propagation and/or rotation 
in interplanetary space. Second, we confirm the previous results 


that the expansion speed is correlated with the radial propaga¬ 
tion speed and most MCs did not expand self-similarly at 1 AU. In 
our statistics, about 62%/17% of MCs underwent a under/over¬ 
expansion at 1 AU and the expansion rate is about 0.6 on average. 
Third, most interestingly, we And that a signiflcant poloidal mo¬ 
tion did exist in some MCs. Three speculations about the cause 
of the poloidal motion are therefore proposed. These findings ad¬ 
vance our understanding of the MC’s properties at 1 AU as well 
as the dynamic evolution of CMEs from the Sun to interplanetary 
space. 


1 Introduction 


Since first identified hY \Burlaaa et ai\ in 1981, magnetic 
clouds (MCs) have been studied extensively in the past 
decades. They are large-scale organized magnetic structures 
in interplanetary space, developed from coronal mass ejec¬ 
tions (CMEs), and play an important role in understanding 
the evolution of CMEs from the Sun to the heliosphere and 
the associated geoeffectiveness. 

The current knowledge of MCs are mostly from in-situ 
one-dimensional observations, and various MC fitting mod¬ 
els have been developed to reconstruct the global picture of 
MCs in two or three dimensions. It is now believed that an 
MC is a loop-like magnetic flux r ope with two ends rooting 
on the Sun [e.g.. 


magnetic tlux r ope witn two ends rooting 
Burlaaa et g^L Il98ll : \Larson et oZ.I . Il997l : 


\ Janvier et al\ . |2013| . The modeling efforts mainly focus on 
two aspects. One is to reconstruct a realistic geometry and 
magnetic field configuration. In past decades, MC fitting 
models have been develo ped from cylindr i cally symmetrica l 


force-free flux ropes [e.g.,[Goldstezrj ,[l98j;|MaroMg4j 


lilirHfflgal . Il988l : Uevvina et al. . 199d : l^irmar and Rust 


198€ 


1991 1 


ropes (e.g.. \Mulliaan and Russell 

. I 2 OOII: \Hu and Sonneruv. 

2 OO 2 I: Hidalao et al.. I20n2a 


Cidetal.. 2002: 

Vandas and RomashetA l200f 

1 and torus-shaped flux ropes 

le.B.. Romashets and Vandas .\200A. Marubashi and Levvind. 

2 OO 7 I: \Hidalao and Nieves-Chinchillah. 2012||. Some com- 


parisons among various MC fitting models could be found 
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in the pap ers by, e.g., \Riley e t al\ 


\Al-Haddad et oZI 


_ iv et aU 

[201l|| a,nd \ Al-Haddad et a/.l [2013| . 

The other important aspect is to understand the ex- 
.nsion and distortion of the cross-section of a MC 


[e.er.. \FarruQia et al\. 

ll99.Sl. I 1995 I: \Marubashi. 

Shimazu and VandaA 

2 OO 2 I: \Berdichevskv et al.. 

Hidalad. l2003l: \Owens et all 20061: \Dasso et al\ 

Demoulin and Dasso. 

2009allbl: Dem,oulin et all 



This is generally an issue about how the velocity is dis¬ 
tributed in an MC. There are lots of observational evidence 
that MCs expand when propa g ating away from the Sun 
[e.g., Klein and Burlaaca. 1 19821: Berdichevsky et al. 1200.1 


I Wano eT'aZri^OOsi r fji^ et al.L l2006l : rCwfeanoTT'^ [ 2 OK 

When the expansion is anisotropic, the initially circular 
cross-section of a MC will deform into a non-circular shape. 
The kinematic treatments and MHD simulations suggested 
that MCs may develop into an ellips e or ‘pancake’ shape 
\IU^ et al\ . 120031 : \Rilev and CroofeCT . I 2 OO 4 I : I Owens et al\ . 
l200fil] . After exploring t he pa rameter space of flux rope 
models, \Demoulin et al\ [2013l| concluded that the aspect 
ratio of the ellipse is around 2, but not too large. On the 
other hand, the fitting results applying an ellipse model 
to the observed MCs suggested that the cross-s ection of 
MCs may be not far from a circle [See Table 1 of \Hidala<i . 
200^. A similar resul t co uld be seen in the pape rs by 


Hu and Sonneru^ |200ll| and lfifn and Sonneruri |2002|] . who 


used the Grad-Shafranov (GS) technique to reconstruct 
MCs in a 2-dimensional plane. This technique does not 
constrain the shape of the MC’s cross-section, but we still 
can find a nearly circular cross-section, particularly for the 
inner core. Thus, it may be acceptable to assume a circular 
cross-section of an MC at 1 AU. 

Furthermore, expansion is only one type of the motion 
of the MC plasma. Imagining a segment of an MC, which 
is approximately a cylindrical flux rope (refer to Fig[T]D), 
one may assume that there are at least three types of the 
global plasma motion: linear propagating motion, Vc = 
(vx ,vy ,vz), expanding motion, Ue and poloidal motion, Vp. 
Vc is the velocity in a rest reference frame, e.g., GSE coordi¬ 
nate system {X, Y, Z), in which spacecraft movement could 
be ignored during the passage of an MC. Ve and Vp are both 
the speeds in a plane perpendicular to the MC’s axis. By 
applying a cylindrical coordinate system, {r,^p,z), sticking 
to the axis of the MC, we have Ve = Vt and Vp = Vp. 

Most of previous studies of fitting locally observed MCs 
simply assumed that MCs propagate radially from the 
Sun, which means that Vc « vxX. However, statistical 
and case studies of the propagation and geoeffectiveness of 



2004 

1 

2011 


agatio n in in t erplanetary s p ace Wm^_et_^ L I 2002 I. 

Ro^Aguez_et_gL 
I 2 OI 4 II . It will obvi¬ 
ously provide a non-radial component of the linear prop¬ 
agation motion. Besides, the orientation of the MC axis 



planetarv space le.g.. R.ust et al.. 20061: Wana et all 2006b 

Yurchvshvri. l2008l: 1 Yurchvshvn et al. 

. 20091: 1 Vourlidas et al. 

2011: Nieve.s-Chinchilla et all 20121: 

Isavnin et all 2014 [. It 


is another source of the non-radial motion. The pictures of 
both deflection and rotation are mainly established on indi¬ 
rect evidence and models. Thus it is interesting to seek any 
signatures of non-radial motion from in-situ data. 


As to the poloidal motion of plasma in MCs, there is 
so far no particular study. But some phenomena hint at 
the possible existence of such a motion. One of them is 
the strong field-aligned s t reams of suprathermal electrons in 
MCs [e.g., \Larson et al ], Il997l[ . Another is the frequently 
observed plasma flows in prominences and coronal loops in 
the solar corona. Prominences and coronal loops may be a 
part of a MC if they are involved in an eruption. If plasma 
poloidal motion did exist in MCs, we will shed new light on 
the dynamic evolution of CMEs. 

The aim of the present work is to investigate the plasma 
motion of MCs from in-situ data with the aid of a flux rope 
model. As the first attempt, we utilize a relatively simple 
and ideal flux rope model, which is cylindrically symmetri¬ 
cal and force-free, but with all the three components of the 
plasma motion taken into account. One could go to the web 
page http://space.ustc.edu.cn/dreams/mc_fitting/ to 
run our model. The details of the model and method we 
employed in this work are described in the next section. 

2 Model and Events 

2.1 Velocity-modified cylindrical force-free flux 
rope model 

2.1.1 Model description 

Our model is developed from the cylindrically symmet¬ 
rical force-free flux rope model that has been widel y used 
in many previous studies [e.g., \Levvina et oZ.I . ll99Cll[ . The 
following coordinate systems are used. One is the GSE co¬ 
ordinate system, {X,Y,Z), in which X is along the Sun- 
Earth (or Sun-spacecraft) line pointing toward the Sun, Z is 
a northward vector perpendicular to the ecliptic plane and 
Y completes the right-handed coordinate system. Another 
is a cylindrical coordinate system in the MC frame, (r, </5, z), 
with z along the axis of the flux rope, as illustrated by Fig¬ 
ure Sometimes to show the pattern of magnetic field and 
velocity, one more Cartesian coordinates in the MC frame, 
{x',y',z'), is used, in which z' is identical with z, x'-axis is 
the projection of the observational path on the plane perpen¬ 
dicular to z', approximately directed toward the Earth and 
y' completes the right-hand coordinate system (see Fig |12b 
and I Ilf 1. 

The magnetic fie ld of a cylindrical flux rope is described 
hv \Lundauisli [l95n| solution in the coordinates (r, ip, z) 


Br 

= 0 

(1) 

Bp 

= HBoJi{ar)(p 

(2) 

B, 

= BoJo{ar)z 

(3) 


in which H = ±1 is the handedness or sign of the helicity, Bo 
is the magnetic field strength at the axis, a is the constant 
force-free factor, and Jo and Ji are the Bessel functions of 
order 0 and 1, respectively. Like common treatment, we set 
the boundary of a flux rope at the first zero point of Jo, 
which leads to a = 2A1/R and R is therefore the radius of 
the flux rope. 

If any kinematic evolution of the flux rope could be ig¬ 
nored, the magnetic field profile along an observational path 
is determined once the following six parameters are known: 
(1) the orientation of the flux rope’s axis, i.e., the z-aixs, 
which is given by the elevation and azimuthal angles, 9 and 
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Figure 1: (a) Schematic picture of an MC at the heliocentric distance of L (adapted from I Wana et aZl [20091] 1. The black 
line indicates the looped axis of the MC with a length of 1. The blue dashed lines suggest the upper and lower limits of 1. 
(b) Illustration of the three types of global motions of an MC. The black, red and blue arrows denote the linear propagating 
motion, expanding motion and poloidal motion, respectively. 


4>, in GSE coordinate system; 6 = 0 is in the ecliptic plane 
and 6 = 0 and 0 = 0 is toward the Sun, (2) the closest 
approach (CA) of the observational path to the flux rope’s 
axis, given by d in units of R; a positive/negative value of d 
means that the observational path is above/below the axis 
in {x', y', z') coordinates, i.e., y' > Q/y' < 0 as shown by the 
examples of the positive d in Fig |12b and 113b . and (3) the 
free parameters in Eq[T}{3l which are H, Bo and R. 

After the possible plasma motion is taken into account, 
five additional parameters should be considered, which are 
Vc = (vx, vy, vz), Ve and Vp. Here we assume that the flux 
rope propagates and expands uniformly and experiences a 
quasi-steady, self-similar expansion (or contraction) during 
the period of interest. Thus, Vc is a constant vector, describ¬ 
ing a global linear motion, that does not change the internal 
magnetic field of the flux rope. 

The parameter Ve is a constant expansion speed of the 
boundary of the cross-section of the flux rope. The self¬ 
similar assumption gives the expansion speed at any radial 
distance, r, away from the flux rope’s axis as 

Vr(x) = XVe (4) 

in which x is the normalized radial distance, equal to r/R. 
As a consequence, the radius of the cross-section of the flux 
rope evolves as follows 


R{t) = R{to) + Ve{t - to) (5) 


in which to is the initial time (or the time the observer first 
encountering the MC flux rope), and the magnetic field as 


Bo{t) 


Bo{to) 


R{to) 

R{t) 


( 6 ) 


Eq[6] follows the magnetic flux conservation. 

It should be noted that, due to the curvature of the whole 
looped flux rope, the propagation velocity of the different 


segment of the flux rope is along the different direction as 
illustrated in Figure [TJd, which may actually cause the ex¬ 
pansion along the axis of the flux rope. That is also why 
the length of the flux rope increases as indicated by Ea ll4l 
below. The axial expansion rate could be on the same order 
of the expansion rate of the flux rope’s c r oss-section ]e.g. , 
Jasso et al\ . l2007l : \DemouUn et al\. I 2 OO 8 I : \Nakwacki et oZ.L 


20111] . which guarantee a roughly self-consistent ev olution 


of a force-free flux rope \Shimazu and Vanda A |2002| . 

Under the assumption of self-similar expansion and the 
conservations in both mass and angular momentum, the 
poloidal speed in the flux rope can be derived (see Appendix 
lATl as 


v^{t,x) = kxfp{x)R(t) 


(7) 


in which fci is a constant and fp{x) is a function dependent 
only on the relative position x. However, the expression of 
fp{x) cannot be specified theoretically. As the first attempt, 
here we tentatively assume fp(x) = 1. It should be noted 
that the point at a; = 0 is a singularity under this assumption 
because U(p(0) has a non-zero value that is not physically 
meaningful. We just ignore this singularity, and as will be 
seen in Sec l4.3[ we prove that this assumption can be treated 
as an acceptable approximation. Then Eq[7|can be rewritten 
as 

vUt,x) = Vp{t)=Vp{to)^^^ ( 8 ) 

where the parameter Vp{t) defines the poloidal speed of the 
plasma at the boundary of the flux rope in the direction of 
0 - 


2.1.2 Parameters 

In summary. Table [T] lists all the 11 free parameters in 
the velocity-modified cylindrical force-free flux rope model. 
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Besides, based on the cylindrical force-free flux rope assump¬ 
tion, we may derive more parameters, which have been also 
summarized in Table [T] It is straightforward to obtain the 
first two derived parameters, tc and 0. The next four pa¬ 
rameters are derived as follows. 

The axial magnetic flux, $ 2 , is given by 

r2Tr pR 

= / / Bzrdrdifi — BqR^ ki (9) 

Jo Jo 

in which ki = ^ xJo{x)dx = 1.35 and xo = 2.41 is the 

Xq 'J'J 

first zero point of Bessel function Jq. The poloidal magnetic 
flux, is given by 

= f f B^drdz = BoRIk 2 (10) 

Jo Jo 

in which I is the length of the loop as denoted in Figure [T]i 
and K 2 Ji(x)dx = 0.416. The helicity, Hm, is given 

by fe.g., \Ber(}er \. |2003ll 

pI p2-k pR 

Hm ~ ^ Jirdrdifdz = BqR^Iks (11) 

Jo Jo Jo 

in which A is a vector potential of B and K 3 = 
25 x[Jq{x) -|- Ji(x)]dx = 0.701. The magnetic energy, 
Em, is given by 

1 pI p2'n pR 

Em = — / / / B^rdrdifdz = BqR^Ika (12) 

2 /^ Jo Jo Jo 

in which m = 6.72 x 10® m H“^. 

We do not know the value of I, but we may assume that it 
is bounded between ttL and 2L, L is the heliocentric distance 
of the flux rope, as shown in Figure [T^. In practice, we let 

The uncertainty, in the length of the flux rope will 

result in the uncertainty in Hm and Em- It should be 
noted that (1) the uncertainty here is not from the fitting 
procedure but the length of the flux rope, and (2) this treat¬ 
ment will underestimate the real length if the leg of an MC 
rather than its leading part was observed. 

Considering the conservation of magnetic flux, we may 
infer from Eg and 1 101 that 

Ron I or R oc L (14) 

and consequently, the helicity is conserved too. Further, 
we may infer Em oc L~^, suggesting that magnetic energy 
continuously decreases when an MC propagates away from 
the Sun. In this work, we calculate its initial value, EmO, as 
listed in Table [U which is the magnetic energy when the MC 
is only one solar radius away from the Sun, and the associ¬ 
ated uncertainty of about ±22% comes from the uncertainty 
in L Be aware that the decay index of the magnetic energy 
of an MC may not be —1. The case study hv \Nakwacki et au 
[^, for example, has shown that the decay index for the 
MC observed by two spacecraft at 1 and 5.4 AU, respec¬ 
tively, on 1998 March is about —0.9, a little bit slower than 
the expectation from our model. If the decay index varies 
from —0.9 to —1.1, the extrapolated EmO suffers an addi¬ 
tional uncertainty of about ^^ 1 %’ "'hich we do not taken 
into account in our model. 

The last derived parameter, Xn, is used to evaluate the 
goodness-of-fit, which will be introduced in the following sec¬ 
tion. 


2.1.3 Evaluation of the goodness-of-fit 

Our fitting procedure is designed to fit the observed mag¬ 
netic field, B°, and velocity, v°, together, which means that 
all the three components of B° and v° are used to constrain 
the model parameters. There are 11 free parameters, among 
which three of them are time-dependent (see Table [T]) and 
can be determined by EqO[6]and[8l Here, we use the first 
contact of the MC, to, as the reference time. The detailed 
fitting process is as follows. First, by given a set of 11 free 
parameters, the imaginarily observational path relative to 
the axis of the flux rope can be derived from the orientation 
of the axis, {6,4>), the closest approach, d, the propagation 
velocity, {vx,vy,vz), and the expansion speed, Ve. Second, 
the coordinates of the observational path are then trans¬ 
formed from the GSE coordinate system to the MC frame, 
{r,ip,z). In the MC frame, the three components of mag¬ 
netic field along the path can be determined by Eq[TU51 in 
which the free parameters Bo, R and H are involved, and 
the velocity can be calculated as Vrr + v,p(p, in which Vr and 
Vif, are given by Eq[4]and[8l Third, we transform the derived 
magnetic field, B"*, and velocity, v"*, in the MC frame back 
to the GSE coordinate system, and evaluate the goodness- 
of-ht by calculating the normalized root mean square (rms) 
of the difference between the modeled and observed values 
of magnetic field and velocity, which is given as 




\ 


1 ^ 

— y 

2N ^ 




± Xvn) 



2 


(15) 


where N is the number of measurements, and Uref is a refer¬ 
ence velocity. 

To find the best-ht, we use the IDL package, MPFIT (refer 
to htt p://purl.com/net/mpfit), to perform least-squares 
fitting \MarkwardA . l2009l : I Morel . Il978|| . The initial value of 
Boijo) is set to be the maximum value of the magnetic field 
strength during the interval of interest, the initial value of 
R(to) is estimated as ^vxAt, in which At is the duration of 
the MC interval and vx is the mean value of the observed 
vx, the initial value of the propagation velocity (vx, vy,vz) 
is the mean value of the observed velocity, the initial value 
of the expansion speed is estimated from the slope of the 
observed radial velocity, and the initial value of the poloidal 
speed is set to be zero. For the free parameter H, we just fix 
its value to 1 or —1 by adding a loop in our fitting procedure. 
The elevation and azimuthal angles, 6 and 4>, a-re two most 
important free parameters. In order to get the best fitting 
result, we test the initial value of 6 every 15° from —90° to 
90°, and do the same thing for (j> from 0° to 360°. Besides, 
we assume that the front and rear boundaries of an observed 
MC define the interval of the flux rope, and then the clos¬ 
est approach, d, could be uniquely determined based on the 
preset velocity and the axis orientation of the flux rope. In 
summary, we try 576 attempts of fitting (i.e., 576 sets of the 
initial values of the free parameters) for an MC, and select 
the case with the smallest value of Xn as the best-fit. 

The reference velocity in Ea ll5l is used to adjust the weight 
of the velocity in evaluating the goodness-of-fit, and is set as 


max(|v°|) -min(|v°|) ^ max(|B°|) - min(|B°|) 
max(|v°|) — Uref max(|B°|) 
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Table 1: Parameters involved in the velocity-modified cylindrical force-free flux rope model’ 


Parameter 


Explanation 


Bo{t) 

m 

H 

e 

</> 

d 

Vx 

Vy 

Vz 

Ve 

Mi) 


tc 

e 

Hm 

E/771O 


X'n- 


Free parameters in the model 
Magnetic field strength at the axis of the flux rope. 

Radius of the cross-section of the flux rope. 

Handedness or sign of helicity, must be 1 (Right) or —1 (Left). 

Elevation angle of the axis of the flux rope in GSE. 

Azimuthal angle of the axis of the flux rope in GSE. 

The closest approach of the observational path to the axis of the flux rope. 

Propagation speed of the flux rope in the direction of X. 

Propagation speed of the flux rope in the direction of Y. 

Propagation speed of the flux rope in the direction of Z. 

Expansion speed of the boundary of the flux rope in the direction of f. 

Poloidal speed at the boundary of the flux rope in the direction of ip. 

Other derived parameters from the model 
The time when the observer arrives at the closest approach. 

Angle between the axis of the flux rope and X-axis. 

Axial magnetic flux of the flux rope. 

Poloidal magnetic flux of the flux rope. 

Magnetic helicity of the flux rope. 

Initial magnetic energy, i.e., the magnetic energy of the flux rope when it was one solar radius away from 
the Sun. 

Normalized root mean square (rms) of the difference between the modeled results and observations. 

See Sec l2.ll for the definition of coordinate systems, derivations of some parameters and other details. 


If there was no reference velocity, XBn and Xvn may not have 
the same weight and cannot be inserted into one formula, 
because the dynamic range of the velocity is much different 
from that of the magnetic field. For example, assuming an 
MG interval during which the magnetic field varies from 5 
to 35 nT and the velocity varies from 300 to 600 km s“^, 
and at a given point B° = 25 nT, S’" = 20 nT, v° = 500 
km s“^ and u"* = 450 km s“^, we can get that the relative 
error between the modeled and observed values in magnetic 
field is ^ = 20% and that in velocity is ^ = 10%. The 
values of the relative errors are different. But if considering 
the range of 30 nT in magnetic field and the range of 300 
km s“^ in velocity, one can find that the relative difference 
between the modeled and observed values in magnetic field is 
the same as that in velocity. Thus the value of relative error 
depends on the dynamic range. To remove this effect, we 
use the reference velocity, which is 250 km s“^ in the above 
case. The corrected relative error of the velocity becomes 
500^250 ~ 20%, the same as that of the magnetic field. After 
this treatment, XBn and Xvn can be fitted into one formula, 
EgfT^ to assess the goodness-of-fit. 

The normalized rms, Xn, has a definite meaning here. It 
measures the average relative error of modeled vectors, B 
and V, with both length and direction taken into account. 
It should be noted that the normalized rms here is dif- 
ferent from those u sed i n some studie s . Fo r example, in 
\Levvina et al\ |l990l| and lLeonmo et oZl [2006| . B° and B"* 
are normalized by the strength of themselves, and that rms 
merely reflects the average de viation of vector directions; 
while m \Marubashi etM\ |2012l| . B° and B"* are normalized 
by the maximum value of observed B during the whole in¬ 
terval of the flux rope, which might imply that a stronger 
MG has a smaller value of rms. 


2.2 Events and model testing 

The MG list compiled by \Levvina et al\ |2006l] is the 
basis of this study, and hereafter we call it Lepping list (see 
http://lepmfi.gsfc.nasa.gov/mfi/mag_cloud_Sl.html 
the list including the fitted parameters are kept being 
updated until 2011 December 13). We use this list not only 
because it is well established but also because we can test 
the procedure of our mod el by comparing our fitting results 
with iLeppfno et~al\ |2006l| results. There are a total of 121 
MCs in the list. The start time, end time and model-derived 
parameters are all given in the list. The goodness-of-fit 
is estimated by Qq, which is 1, 2 or 3 for good, fair and 
poor, respectively. Our study only considers the MCs in 
the first two categories. It is noticed that there are large 
data gaps in the published Wind data during 2000 July 
15 “ 16. Their event No.45 and 46 are in that period, and 
are therefore removed from our study. Besides, th e event 
No.85 in their list was studied bv lDasso et al\ [2009l| . and is 
believed to consist of two MCs. Thus we discard this event 
too. Finally, a total of 72 events are available for our study 
as listed in Table EMI in Appendix m 

To test our procedure, we switch off the velocity option 
in our model, and compare the fitted parameters with those 
given in Lepping list. When fitting an observed MG with 
our model, we use the time-resolution of the Wind data as 
the same as that indicated in the list, which is 15 minutes for 
some events and 30 minutes for others. The front and rear 
boundaries of the MCs given in Lepping list are used and 
fix ed in our model. It is different from the procedure taken 
hv \Levvino et al\ [2006| . in which the preset boundaries may 
be modified by their model. The difference between their fit¬ 
ted and observed boundaries is evaluated by the parameters. 
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Figure 2: Histograms showing the differences of the values of fitting parameters between the Lepping’s model and our 
model without velocity measurements taken into account. From top to bottom and left to right, there are the relative 
difference of the magnetic field strength (Bo) at the MC’s axis between the two models, the relative difference of the radius 
(R), the difference of the closest approach |ci|, the acute angle between the MC’s axes derived from the two models, and 
the difference of the handedness. 








ck an d asf, given in their list (see Eq.7 and 9 m \Levvina et al\ 
l200fil for their definitions). Larger values of ck and asf mean 
larger difference. Thus to make a mostly reliable compari¬ 
son, we further tentatively exclude events whose ck or asf is 
larger than 10%, that results in a sample of 23 events. It 
should be noted that those excluded events will be included 
again in the next section though they are not used for testing 
our procedure here. 

The parameters characterizing the configuration and 
strength of an MC are compared, which are H, Bo, R, d, 
and orientation {9 and 0) of an MC. Figure [2] shows the 
diffe rences of the values between the parameters from our 
and lLeppmo et oil [200fil] procedures. It is found that both 
of them suggest the same sign of helicity, H, for each MC 
(Fig.[2]3). For Bo and R, we calculate the relative difference, 
where / is from our results and /l from Lepping list. 
Our procedure gets very similar results with Lepping list 
(Figl2^ and[2j3)- There are only 3 events, in which the rela¬ 
tive difference in Bo and R is outside of ±20%. For d, we do 
not use the above equation to calculate the relative differ¬ 
ence, because d is a ratio of the closest approach to the radius 
of the flux rope and its value could be zero which may cause 
the value of the relative difference to be extremely large. 
Instead, we directly use the difference of absolute values of 
them, |d| — \dL\- A good agreement also can be obtained 
between the two procedures (Fig[2j:). The two parameters 
9 and 4>y which define the orientation of a flux rope’s axis, 
are compared together. We calculate the angle between the 
orientation derived by our procedure and that given in Lep¬ 
ping list. It is found that the angle is smaller than 60°, and 
for most (19 out of 23) events, the angle is less than 30°. 
Overall, the comparison shows that our procedure can al - 
most reproduce the results derived '\yv \Levmna et aZl |2006l| . 
conhrming the validity of our procedure. 

Now we switch on the velocity option, and apply the model 
to the 72 high-quality events in Lepping list. The time- 
resolution of the data input into our model is set to 10 min¬ 
utes. All the parameters derived from our htting procedure 
of these events have been listed in Table [2] [4] in Appendix m 
for reference. In the following two sections, we will evaluate 
the effects of velocity on the goodness-of-fit, and show the 
statistical properties of plasma motion. 

3 Effects of velocity on the fitting results 

The goodness-of-fit is evaluated by the average relative 
error, Xn, given by Eg 1 151 For all the 72 events, the distri¬ 
bution of Xn has been displayed in Figure [Hfi- It is found 
that the values of Xn are all less than 60%, and on average, 
Xn is about 28%. Further, the comparison between the cases 
of velocity-on and velocity-off shows that considering veloc¬ 
ity may get the value of XBn larger or smaller (see Figl^fi). 
As examples, Figure IHi and [IJs show two cases; for one of 
them the fitting to the magnetic field gets worse, and for 
the other, the fitting looks better. In a typical MC measure¬ 
ments, a declining velocity means the expansion of the MC 
and as a consequence, the magnetic field strength of the MC 
should decrease as well. The former case does not follow the 
picture. The strongest magnetic field appeared in the rear 
part of the MC, but the velocity profile indicates a clear ex¬ 
pansion, suggesting that the strongest magnetic field should 
appeared in the front part of the MC. The fitting without ve- 




Xsn XsnO 

Figure 3: The distributions of normalized rms, which show 
the goodness-of-fit. (a) the normalized rms (Xn) calculated 
from both magnetic field and velocity, and (b) the difference 
of the normalized rms in magnetic field between the model 
results with velocity-on (xBn) and velocity-off {xBno)- 

locity being taken into account gives the blue dashed curves, 
which match the observed magnetic field much better than 
red solid curves given by our velocity-modified model. The 
latter case is a typically expanding MCs, and therefore the 
magnetic field is fitted better. 

Except those velocity-related parameters (which will be 
presented in the next section), the parameters of MCs de¬ 
rived from our velocity-modified model are summarized in 
Figure [S] For the magnetic field strength at the MC axis, 
most events fall into the range from 10 to 30 nT with an 
median value of about 16 nT. For the radius, all the MCs 
are less than 0.25 AU, and its median value is about 0.09 
AU. Note, the values of the above two parameters are all 
adopted at the time, tc, when the observer arrives at the 
closest approach to the MC axis. The angle between the 
MC axis and the Sun-spacecraft line, 0, could be any value 
from 0° to 90°. The most probable angle is within 45° — 75°, 
and the median angle is about 53°. It suggests that observed 
MCs are more likely to transversely cross over the observer. 
The elevation angle of the MC axis tends to be small, with 
a median value of about 15°, implying low tilt angles of flux 
ropes when they erupted from the Sun. 

Further by using the total length of the flux rope, I, given 
by Ea ll3l we estimate that the total magnetic flux, which in¬ 
cludes axial flux and poloidal flux, is on the order of 10^^ Mx 
with the median value of about 4.1 x 10^^ Mx. The distri¬ 
butions of the axial flux and poloidal flux are also indicated 
by dotted and dashed lines (Fig[S]B). Their median values 
are 0.4 x 10^^ and 3.6 x 10^^ Mx, respectively. The estima- 
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Figure 4: Two cases showing the comparison of the fitting results of our model with velocity-on (the red lines) and off 
(the blue dashed lines). Case (a) shows that the fitting results become worse when the velocity is taken into account, 
and Case (b) shows an opposite situation. In both cases, from the top to bottom, the panels are total magnetic field 
strength, elevation and azimuthal angles of magnetic field vector and three components of bulk velocity of solar wind 
plasma, respectively. The two vertical lines mark the front and rear boundaries of the MCs. 


tion of the magnet i c flux of MCs were mad e before by, e.g., 
IfJosso et oL 2005|[j_ Dassoet al\ [20071] and \Nakviacki et al\ 
[2008l| . In iDasso et al\ [2005l| . for example, eight well-defined 
MCs were investigated and it was found that the axial flux 
is around 0.4 x 10^^ Mx, highly consistent with the result 
obtained here. Our results are also roughly in agreement 
with previous studies about the magnetic flux of MCs an d 
reconnection flux of solar eruptions by, e.s... \Om et al ]1^. 

The helicity is shown in Figure [3' and [SJi. The number 
of right-handed MCs is almost equal to the number of left- 
handed MCs, and the absolute value of the helicity is about 
2.05 X 10**^ Mx^ o n average. For the 8 events studied by 
iDasso et al\ [2005| . the helicity per unit length is about 1 x 
lO'^^ Mx^/AU. According to Eq[T3l the helicity in their study 
is about 2.57 x 10“*^ Mx^, consistent with the average value 
we obtained here. Moreover, the initial magnetic energy is 
found to be about 2.84 x 10^^ erg, nearly one order higher 
than the kinetic en ergy of a typical CME , which is about 
10^® — 10^° erg \e.e:.. \Vourlidas et aLLl20iri| . even considering 
the uncertainty in the decay index when we extrapolate the 
in itial magnetic energy. 

\Nakwacki et al\ [20081] compared the static and expansion 
models and found that the modeled values of the MC pa¬ 


rameters are not changed too much. Here we consider not 
only the expansion but also other types of motion. Will the 
values of parameters change more significantly? By compar¬ 
ing the values of the fitting parameters obtained from the 
velocity-modified and non-velocity models as shown in Fig¬ 
ure [6| we find that the values are more or less changed, but 
not too significantly except the orientation. The magnetic 
field strength is almost unaffected. For the radius, about 
13% (9 out of 72) of the events get a smaller value after 
velocity is taken into account, and about 10% of the events 
get a larger value. For the closest approach, the velocity- 
modified model believes that the observer should be slightly 
farther away from the MC axis for 8 events and closer for 
only 4 events. The largest difference appears in the orien¬ 
tation. There are about 40% of the events, for which the 
orientation is changed by more than 15°, and particularly, 
there are 6 events with the difference in orientation larger 
than 45°. Besides, for two cases, the handedness is also 
changed. 
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Figure 5: The distributions of the values of the fitting parameters from the velocity-modified model for 72 MC events. 
From the top to bottom and left to right, the panels show the magnetic field at the MC’s axis (Bo), the radius the MC’s 
cross-section (R), the acute angle between the axis and the Sun-spacecraft line (0), the absolute value of the elevation 
angle of the axis (|0|), the total magnetic flux ($2 + 'I’,,^), the unsigned helicity (|7/„,|), the initial magnetic energy (Emo) 
and the handedness (H). The median values of these parameters are indicated by the arrows. 
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Figure 6: Histograms showing the differences of the values of fitting parameters between our model with and without 
velocity measurements taken into account. The arrangement is as the same as that in Fig. [2] 
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4 Statistical properties of plasma motion 

4.1 Linear propagating motion 

One important implication to consider the velocity in fit¬ 
ting procedure is to reveal the properties of the plasma mo¬ 
tion of MCs. First, we investigate the axial component, 
Vaxia, and perpendicular component, Vperp, of the linear 
propagating motion of an MC. The former is a velocity along 
the MC’s axis, i.e., in the ^-direction of the MC frame, and 
the latter is a velocity perpendicular to the radial direction, 
i.e., the Sun-spacecraft line, Vperp = Vvy + ^1- 

The ratio of the axial velocity to the radial velocity, 
—Vaxis jvx, as a function of the orientation of the MC’s axis 
is shown in Figure [Tli. A very nice correlation between them 
could be found. When the axis is almost aligned with Sun- 
spacecraft line, the absolute value of Vaxis is almost equal 
to that of vx, and when the axis becomes more and more 
perpendicular to the Sun-spacecraft line, Vaxis approaches 
zero. It well follows the cosine function as indicated by the 
solid line. The small deviation away from the cosine function 
is no more than O.Ohnx, which corresponds to a very small 
speed on the order of 10 km s“^. This result suggests that 
the apparent axial velocity is mostly a consequence of the 
propagation of the MC though an insignificantly pure-axial 
flow might exist inside an MC. 

The distribution of Vperp is given in Figure [S] Panel (a) 
shows that except 6 events, all the other events have a per¬ 
pendicular velocity less than 60 km s“^, and Panel (b) sug¬ 
gests that in 61 (about 85% of) events the perpendicular ve¬ 
locity is no more than 10% of radial velocity. A noteworthy 
thing is that some events have a signihcantly perpendicular 
velocity with respect to the radial velocity. The MC shown 
in Figure|4]D is an example, which occurred on 1998 Septem¬ 
ber 25. Our model infers that the orientation of the MC 
axis is 0 = 60° and tj) = 196°, and the plasma motion of the 
MC is the combination of the linear motion at —624, 95 and 
50 km s“^ in X, Y and Z directions and the expansion and 
poloidal motion at 90 and —23 km s“^, respectively. 

Considering the expansion of the whole looped structure 
as shown in Figure[Tj3 and that there is no signihcantly pure- 
axial flow inside an MC, one may hnd that a signihcantly 
perpendicular velocity may be present if not the front part 
(or the apex) but the hank (or the leg) of the MC was de¬ 
tected locally. Which part of a looped MC is detected by 
spacecraft could be roughly inferred from the orientation of 
the MC axis. We think that the apex of an MC is encoun¬ 
tered if the inferred axis of the MC is almost perpendicular 
to the Sun-spacecraft line, and the leg is encountered if the 
axis is almost parallel to the Sun-spacecraft line. Thus, if the 
expansion of the whole looped structure was the reason of 
the presence of the perpendicular velocity, we may expect a 
very small perpendicular velocity for those apex-encountered 
events. Figure[7l3 shows the scatter plot of —Vparp/vx versus 
0. Obviously, there is no dependence of the perpendicular 
velocity on the axis orientation. For those events with 0 
close to 90°, the perpendicular velocity could be larger than 
10% of the radial velocity. Thus, there must be some other 
reasons for the signihcantly perpendicular motion. 

We think that the presence of a signihcantly perpendicu¬ 
lar velocity might be an in-situ evidence of de hected propa- 
gation of a CME in interplanetary space [e.g.. I Wana et ali . 


|2004 |2014 \Luaa^. 120101: \lsavnin et al. I EH. The positive 
value of vy of the 1998 September 25 event suggests an east¬ 
ward dehection in the eclipt ic plane. This is in agreement 
with the proposed picture bv l Wana et al\ |2004l| that a CME 
faster than background solar wind will be dehected toward 
the east. If a CME propagated in interplanetary space with 
a perpendicular velocity at a tenth of the radial velocity, i.e., 
dehection angle is then given by 




'^perp 

vxL 


dL 



(17) 


Assuming Li = 1 AU and Lq less than 5 solar radii, we can 
estimate that the dehection angle is more than 20°. This 
is quite consistent with our recent case study of a CME, 
which was found to be de hected by mo r e tha n 20° on its way 
from the corona to 1 AU [WanQ et al\ . l2014|| . Alternatively, 
the perpendicular velocity might also be the result of the 
rotation of the whole structure of a CME with respect to the 
radial direction in interplanetary s pace, just as the rotation 
in the middle and outer corona \e .s..\Yu rchvshvn et aLLl2009l : 
I Vourlidas et al. I EH; I Jsaumn et al.l I2OI4I] . 


4.2 Expanding motion 

The distribution of the expansion speed is shown in Fig¬ 
ure [3i. It is noteworthy that a signihcant fraction (about 
26%) of the events experienced a contraction process with 
a median value of about 12 km s“^. We check the large 
contraction events having Ve < —20 km s“^, and hnd un¬ 
surprisingly that they were all caused by the overtaking of 
faster solar wind stream, as illustrated by the example in Fig¬ 
ure nni In that event, the solar wind speed after the trailing 
edge of the MC is much larger than that before the leading 
edge, which cause the magnetic held strength increase with 
time, and reach the maximum near the rear boundary. Un¬ 
der this circumstance, the MC cannot expand freely, but will 
be compressed by ambient solar wind. Our model suggests 
that the contraction speed of the MC is about 27 km s“^. 

For the rest of the events suggested to be expanding at 1 
AU, the median speed is about 21 km s“^. The expansion 
substantially depends on the balance between the internal 
and external forces. The internal force is partially char¬ 
acterized by magnetic held strength. Figure [^3 shows the 
correlation between Ve and Bq. The black dashed line is 
the linear htting to all the data points and the blue solid 
line is the linear htting to the data points of Ue > 0 (blue 
dots). A weak correlation could be found, but it is not sig¬ 
nihcantly different between all data points and data points 
of Ve > 0. The correlation coefhcient is not so high, because 
the external condition is not considered. With the increasing 
distance away from the Sun, both the magnetic and thermal 
pressures in the solar wind usually decrease. It will conse¬ 
quently cause the unbalance between the internal and ex¬ 
ternal forces, which drives a CME/MC expanding. Thus we 
may expect that the expansion speed should be correlated 
with propagation speed. As shown in Figure [HJ;, there does 
exist a stronger correlation between Ve and vx for expansion 
events. The correlation coefhcient is a bout 0.55. It is consis¬ 
tent w ith the previous result by, e.g., \Demoulin and Dassd 
|2009a| that rapid decrease of the total solar wind pressure 
with heliocentric distance is the main driver of the MC ex- 
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Figure 7: (a) The axial velocity and (b) the perpendicular velocity as a function of the angle between the MC’s axis and 
the Sun-spacecraft line. Both the velocities are normalized by the value of radial propagation velocity. A small value of 0 
means an encounter of the leg of an MC and a large value of 0 means an encounter of the apex of an MC. 
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Figure 8: Distributions of (a) the propagation velocity perpendicular to the radial direction, i.e., the Sun-spacecraft line 
and (b) its value relative to the radial propagation velocity. Median value is indicated by the arrow in each panel. 


pansion. Similar results could a lso be found for the events in 
the inner and outer heliosphere \Gulisano et aZ.I . l2010l . l2012| . 

On the other hand, the MC expansion could be classi¬ 
fied as self-similar expansion, overexpansion and underex¬ 
pansion. Note, the term ‘self-similar’ here only refers to that 
the radius, R, of an MC evolves proportionally to the dis¬ 
tance L, which is a subset of that defined in Sec.2.1 where 
‘self-similar’ means that not only the size but also the in¬ 
ternal plasma parameters of the MC evolve self-similarly. 
Imaging data have suggested that most CME s undergo a self¬ 
simil ar expansion in the outer corona [e.g., \Schwenn et oZ.I . 
|200!tI] , but whether or not they maintain the self-similar 
expansion in interpla netary space? To m easure the ex¬ 
pansion rate of MCs, \ Gvlisano et al ] [2nin| used a dimen¬ 
sionless quantity which Aux is the dif¬ 

ference of the measured solar wind speed along the Sun- 
spacecraft line between the front and rear boundaries of 
an MC. They found that the value of is about 0.7 af¬ 
ter analyzing all the MCs observed by Helios spacecraft. 
Considering Aux is a proxy of the expansion speed of an 
MC and vxAt approximates the size of the MC, we may 
infer that C, ~ and ^ < 1 means an underexpan¬ 

sion. Thus, the parameter, C,, has the same physical mean¬ 
ing of the power index, n, appearing in the power law of 
the heliocentric distance dependence of the CME/MC size 


e.g.. \Bo thmer a nd Schwenri . Il998l : \Demoulin et al\. l2008l : 


\Savani et al. . 2009| . na mely S = cL^, where S is th e MC 


size and c is a consta nt. \Bothrner an d Schwenri [l998l] found 


n « 0.78 a,nd \Demoulin et al\ [2008l| gave n « 0.8. Here we 
use model derived parameters to investigate the expansion 
rate again. Since all the MCs in this study located at 1 AU, 
we have L = 1 AU. Figureshows the plot of —v^jvx ver¬ 
sus i?/AU. The self-similar expansion is given by the dashed 
line. The data points above the line suggest an overexpan¬ 
sion, and those below the line a underexpansion. By consid¬ 
ering a 20% uncertainty as indicated by the two dotted lines, 
we find that only 21% of the events underwent a nearly self¬ 
similar expansion, and 62%/17% of the events have and ex¬ 
pansion rate lower/larger than 0.8/1.2. By using a function 
of — CtXu these data points as indicated 

by the solid line, we get / « 0.6 on average, generally con¬ 
sisten t with that obtained by , e.g., Botivm&r^^iidSdwiem^ 
[l998l| , \Demoulin et al\ [2008f| and \ GulisanoTr^ill ' \201(j^ . if 
the uncertainty is considered. A possible reason of why the 
expansion rate is significantly below unity is that the MCs 


are perturbed by ambient solar wind and /or other tran¬ 
sient s as shown in observational statistics \Gulisano et aZ.I . 
l201Cl|| . which may cause the external pressure surrounding 
the MC decreasing with distance more slowly than usua l. 
The numerical sim ulations by, e.g.. lAiono et all [20061 . 120071] . 
\Luoaz et al\ [2018| also suggested that CME-CME interac- 
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Figure 9: (a) Distribution of the expansion speed. The two arrows indicate the mean values for negative and positive 
expansion speeds, respectively, (b) Scatter plot showing the correlation between the expansion speed and magnetic field 
strength. Blue dots indicate the expansion events and black dots the contraction events. The black dashed and blue 
solid lines are the linear fitting to the all data points and the blue data points, respectively, (c) Scatter plot showing the 
correlation between the expansion speed and the radial propagation speed. The line gives the linear fitting result, (d) 
Scatter plot showing the correlation between the expansion speed in units of the radial propagation speed and the MC’s 
radius in units of AU. The solid line gives the fitting result by a function oi y — ax as marked in the upper-left corner of 
the panel. The dashed line indicates the self-similar expansion, i.e., R evolving proportionally to L, and the two dotted 
lines give the 20% uncertainty. For the panel (c) and (d), only expansion events are included. 


tion may affect the expansion rate of the preceding CME. 

It should be noted that in our model the expansion speed 
is derived based on the measured solar wind velocity along 
the Sun-spacecraft line and therefore the expansion speed 
here is more likely to reflect the radial expansion rather 
than lateral expansion. Lots of studies have shown that 
the cross-section of a CME will more or less di storted from 
the circular shape to a pancake shape [e . g.. \Rilev et aH 
I 2 OO. 1 I : \R.ilev and Crook^ . l2004l : I Owens et o/.l . l200fil] . suggest- 
ing that the lateral expansion is probably faster than the 
radial expansion. Thus, the lateral expansion may probably 
still be self-similar as that in the outer corona though the 
radial expansion is not. 

One may notice that our model implies that R oc L, mean¬ 
ing a self-similar expansion (see Sec l2.1.2'1) . It is inconsistent 
with the result here. This inconsistency may come from vari¬ 
ous sources, e.g., the assumption of the magnetic flux and/or 
helicity conservation, the assumption of L oc Z, and the as¬ 
sumption of the uniformly straight cylindrical geometry. It 
will consequently affect the values of some derived parame¬ 
ters, e.g., Hm and Em- Although such an inconsistency 
exists, we think that the values of these derived parameters 


still could be treated as a first-order approximation. 

4.3 Poloidal motion 

The distribution of poloidal speed is shown in Figure fTTI 
The median value is about 10 km s“^. By comparing the 
poloidal speed with other parameters, we cannot find any 
significant correlation among them (all the correlation co¬ 
efficients are no more than 0.40). Obviously, the poloidal 
speed, if any, is less significant than the expansion speed on 
average. This fact may cause the observational signature of 
the poloidal motion in an MC unclear. To most clearly show 
the poloidal motion, we choose the events without significant 
expansion and convert the velocity into the Cartesian frame, 
{x',y',z'), of the MC (see Sec l2.1.l1 for the definition of the 
coordinates). A nice example is shown in Figure [121 which 
was observed on 2009 October 12. The two vertical fines 
mark the beginning and the end of the MC, and the fitting 
results are given by the red curves in Figure fT^ . According 
to the modeled parameters, the path of Wind spacecraft in 
the MC frame is shown in Figure [T^ . Along the path from 
the beginning to the end of the MC, the rotation of the mag- 
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Figure 10: An example of contraction events. The arrange¬ 
ment is as the same as that in Fig. 01 


netic field vector is well presented in both the (x',y') and 
{y',z') planes (see Fig ll2b and I12b l. The color-coded dots 
are observations and the color-coded lines are fitting curves. 
More interestingly, the poloidal motion is evident in the MC 
frame (Fig[T2ji). In this event, the expansion speed is almost 
zero and the poloidal speed is about —31 km s~^ at time of 
tc- Figure [13] shows another event on 2003 March 20, in 
which the modeled poloidal speed is 58 km s“^, more signif¬ 
icant than the modeled expansion speed of about 19 km s“^, 
and therefore the poloidal motion can be also recognized in 
observations. 

Based on the model, we could expect that the data points 
of velocity will form an arc which is symmetric about the 
axis of Vyi = 0 if the poloidal motion was significant and 
dominant in an MC. The above two cases just show the pat¬ 
tern. However, the arc’s symmetrical axis will change to 
v^i = 0 if the expansion was significant and dominant. Fig¬ 
ure [14^ roughly shows the case, which was observed during 
1998 November 8-10 and the values of and Vp are 69 
and —8 km s“^, respectively. If both expansion and poloidal 
motion are significant, the ‘symmetric’ axis will rotate, and 
the velocity data points may deviate slightly from a sym¬ 
metric distribution. The 2001 April 4-5 event shown in 
Figure mb is an example to show the change. 

It should be noted that fp{x) in the equation of poloidal 
speed (EqjTjl is assumed to be unity. Whether or not is this 
assumption reasonable? We check it by investigating the cor- 
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Figure 11: Distribution of the poloidal speed. 

relation between fp{x) = and x as shown in Figure [T51 
Vtp is derived from velocity measurements in the MC frame, 
and R are obtained from model results. ki is an unknown 
constant, that may change from one event to another. We 
determine the value of fci by making the value of fp{x) of 
most data points around unity. All the data points of the 
72 events are plotted together to show the statistical trend 
IFig llSh '). The red diamonds in the hgure are the mean 
values of the data points in the bins with a horizontal size 
of 0.2, and the error bars indicate the standard deviations. 
The blue line gives fp(x) = 1. It is clear that data points are 
generally distributed around the blue line no matter which 
value X is, though a large scattering is evident, suggesting 
that fp{x) = 1 is an appropriate assumption from the view 
of statistics. 

However, the form of fp{x) is case dependent. Is the as¬ 
sumption of fp{x) = 1 generally appropriate for an individ¬ 
ual MC? It is further checked in Figure I15b415h . in which 
the dependence of fp{x) on x for 7 selected MCs are pre¬ 
sented. These 7 MCs are selected for investigation because 
(1) the derived poloidal speed is signihcant, larger than 10 
km s“^, and (2) the closest approach is smaller than 0.57? 
so that there is a complete scan of x. For different events, 
the pattern is indeed different. Figure fT5b and I15h suggest 
a general increase of fp{x) with increasing x, Figure [T3h and 
EE show that fp{x) has a unimodal distribution, and Fig¬ 
ure EEiEE and Elk show a reversed unimodal distribution 
of fp{x). If considering the errors, we find that fp{x) is al¬ 
most independent on x for all the cases except the last one 
(FigEEi). Even for the last case, fp(x) is almost invariant 
for X > 0.4. Thus, based on the above analysis, we may 
treat the simplest assumption, fp(x) = 1, as an acceptable 
approximation. 

5 conclusions and discussion 

In this paper, we present a velocity-modified cylindrical 
force-free flux rope model in details. Both observations of 
in-situ magnetic field and plasma velocity are taken into ac¬ 
count by our model to derive the geometrical and kinematic 
parameters of MCs, which have been summarized in TableEl 
The validity of our fitting procedure and the effect of velocity 
on the fitting results are checked through the Lepping list. 
It is found that the values of the modeled radius and orien- 
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Figure 12: An MC observed on October 12, 2009. (a) The fitting results of the velocity-modified model. The arrangement 
of the panels is as the same as that in Fig. [T] (b) and (c) The magnetic field components in the [x', y') and (y', z') plane of 
the MC frame, respectively. The data points are color-coded, indicating the time from the beginning (black, corresponding 
to the hrst vertical line in the plot (a)) of the MC to the end (red, corresponding to the second vertical line in the plot 
(b)). The color-coded lines are the fitting results, (d) The velocity in the (x',y') plane, (e) The derived cross-section of 
the MC (the circle) and the observational path (the arrow). See main text for more details. 


tation of an MC and the closest approach are more likely to 
be changed significantly if the velocity is considered. In our 
sample, the radius changes its value by more than 20% in 
22% of the cases, orientation changes by more than 30° in 
15%, and the closest approach changes by more than 20% in 
17%. In a few cases, the handedness may also be changed. 
We then obtain the statistical properties of MCs, including 
the magnetic field strength, radius, orientation, magnetic 
flux, helicity and initial magnetic energy, which have been 
summarized in Figure [S] Furthermore, some findings about 
the plasma motion of MCs are obtained. 

1. The linear propagation velocity may not be along the 
radial direction. The value of the non-radial component 
of the propagation velocity could be more than 10% of 
radial speed in some cases, which constitutes the direct 
evidence of the deflected propagation and/or rotation 
of a CME in interplanetary space. 

2. As previous studies have shown (see Sec l4.2ll . the ex¬ 
pansion speed is correlated with the radial propagation 
speed with a coefficient of about 0.55, and most MCs did 
not undergo a self-similar expansion at 1 AU in radial 
direction, i.e., the radius is not evolving proportionally 
to the heliocentric distance. In our statistics, 62%/17% 
of MCs underexpanded/overexpanded with an expan¬ 
sion rate < 0.8/ > 1.2, and on average the expansion 
rate is about 0.6. 


3. The poloidal motion is not as significant as expanding 
motion generally, but does exist in some cases. Its speed 
is on the order of 10 km s“^ at 1 AU. 

The last point is of particular interest. Considering the 
possible presence of small pure-axial velocity (see Sec ld.lll . 
there may exist helical plasma motion following the helical 
magnetic field lines of an MC, which is worth to be inves¬ 
tigated further. No matter how the plasma elements move 
inside the MC, the presence of the poloidal motion means 
that MCs may carry non-zero angular momentum. We think 
that there are at least three possible causes of the angular 
momentum. 

The hrst is that the angular momentum is generated lo¬ 
cally through the interaction with ambient solar wind. If 
there is the velocity difference between the solar wind and 
an MC, the solar wind plasma would stream around the MC 
body. The viscosity or some other processes may cause the 
poloidal motion inside the MC. If this is true, we would ex¬ 
pect that the poloidal speed, Vp, is stronger at the periphery 
of the hux rope and therefore the absolute value of the de¬ 
rived Vp is positively correlated with the closest approach, 
d. Our statistical study suggests that there is a weak corre¬ 
lation between them as shown in Figure [TBh: the correlation 
coefficient is about 0.4. This result implies that the local 
solar wind interaction perhaps is probably a factor for the 
poloidal motion but not the only one. 
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Figure 13: Similar to Fig. [12] but for an MC observed on March 20, 2003. 


The second possible cause is that the angular momentum 
is generated internally all the time of the CME from its 
early phase in the corona to propagation phase in interplan¬ 
etary space. It is well known that the m agnetic energy of a 
CME will decrease when it expands [e.g. J JfMmor and Rusu . 


llOOOl : I Wana et al\ . l2009l : \Nak'wacki et alV 20 111) . The de- 


creased magnetic energy may go into thermal energy and 
kinetic energy through some mechanisms. It means that the 
poloidal motion could be somehow generated inside the MC. 
If this is true, the decreased magnetic energy must partially 
go into the rotational kinetic energy, and we would expect 
that the poloidal speed, Vp, be more or less dependent on the 
radial propagation speed, vx, because generally the faster 
the propagation speed is the larger is the rate of the mag¬ 
netic energy decrease. As we can see in Figure fTOb . there is 
also a weak correlation with a coefficient of 0.4. It suggests 
that the internal process perhaps is another factor for the 
poloidal motion. The nu merical simulation s of the evolution 
of magnetic flux rope by I Wei et al\ |l99ll] had shown that 
the expansion of a flux rope may cause the significantly az¬ 
imuthal velocity, i.e., the poloidal speed in our study, even 
if the azimuthal velocity was zero initially. 


The third one may be that the angular momentum is gen¬ 
erated initially at the eruption of the CME in the corona, 
and carried all the way to 1 AU. The kinking/unkinking be¬ 
havior and rotation of the whole erupting magnetic structure 
are often observed in co r onagraphs [e.g., 120081 : 

I Yurchvshvn et aTl l2009l : I Vourlidas et al\ . 20nir Such pro¬ 
cesses might introduce poloidal motion in the body of the 
CME. However, if the angular momentum is conserved as 
that we treated in the model (see Eq0, we may estimate 
that the poloidal speed will be about 200 times of that at 1 
AU, or about 2000 km s“^ near the Sun if the poloidal speed 


is about 10 km s“^ at 1 AU. It is interesting to check the 
imaging data of CMEs to search such a fast poloidal motion 
though no such a phenomenon was reported so far. 

No matter which of the above speculation(s) is right, the 
evidence of poloidal motion inside MCs does shed a new light 
on the dynamic evolution of CMEs from its birth to inter¬ 
planetary space. Some outstanding questions are open. How 
is the angular momentum in the solar wind transfered into 
an MC if the first speculation is right? How is the magnetic 
energy of an MC converted into the rotational kinetic energy 
if the second speculation is right? How does the poloidal mo¬ 
tion depend on the heliocentric distance? Whether or how 
does the poloidal motion modify the macro/micro-scopic 
properties of an MC? Except further studies based on the 
current observations, the upcoming missions, Solar Orbiter 
and Solar Probe+ will provide great opportunities to ad¬ 
vance our understanding of these questions. 
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A Derivation of poloidal speed i;^ 

Under the self-similar assumption, we may assume the 
self-similarity variable x = r/R{t), and the self-similarity so¬ 
lution p{t,r) = po{t)fp{x) and v^{t,r) = vo{t)fp{x), which 
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Figure 14: The scatter plots of the velocity components in the {x',y') plane for (a) the MC observed during November 
8-10, 1998 and (b) the MC during April 4-5, 2001. The arrangement is as the same as that in Fig. 1121 1. 


mean that the shape of the spatial distribution of p and Vy, 
along the radial distance in the flux rope will not change 
with time. Further we have the self-similarity solution of 
Vr{t,r) = Ve{t)x = in which Ue(t) = ^ is the expan¬ 

sion speed and R(t) is the radius of the flux rope. 

The mass conservation requires 


/ 


prdrd(pdz = 27rl 
1 


f 


prdr = 2nlR / pxdx 


/■ 


= 2tvIR^ po{t) [ fp{x)xdx = M 

Jo 


(18) 


Po{t) = 


^-R-^r^ 


27rci 


(19) 


po{t) = 


-R- 


( 20 ) 


/ 


prViprdrdifidz = 2 ttIR I pv^pX dx 


f- 


Vo{t) = 




27rc3 


C 3 M 


( 22 ) 


vo(t) = kiR 


(23) 


where fci = pipf- is also a constant, and therefore expression 
of Vp writes kifp{x)R~^, as written in Eq[71 


B Additional tables 


where M is the total mass of the flux rope, and the above 
integral is a constant, say ci. Then we have 


The conservation of fluxes implies that I oc R fEa ll4ll or 
I = C 2 R, and therefore the above equation becomes 

M 


27rciC2 

Similarly, if there is no force in (p direction, the conserva¬ 
tion of angular momentum requires 


2'kIR^ pQ{t)vo[t) f fp{x)fp{x)x‘^dx = La ( 21 ) 

Jo 

where La is the total angular momentum, and the integral 
in the above equation is also a constant, say C 3 . Then we 
get 
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Figure 15: Scatter plot of measurements during the MC events showing the validity of the assumption of fp{x) = 1 in 
EqEl Panel (a) contains all the data points of the 72 MC events, and each of the other panels display the data points in 
a selected MC. See the main text for more details. 




Figure 16: (a) Scatter plot showing the correlation between the poloidal speed and the closest approach and (b) that 
between the poloidal speed and radial propagation speed. The lines give the linear fitting results. 
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Table 2: Parameters of 72 MCs derived by the velocity-modified cylindrical force-free flux rope model 
MC Interval I Modeled Parameters 


No. 

(1) 

^0 

(2) 

At 

(3) 

Bo 

(4) 

R 

(5) 

H 

(6) 

e 

(7) 

00 

d 

(9) 

vx 

(10) 

Vy 

(11) 

Vz 

(12) 

We 

(13) 

Vp 

(14) 

Ate 

(15) 

e 

(16) 

(17) 

(18) 

Hm 

(19) 

EmO 

(20) 

Xn 

(21) 

1 

1995/02/08 05:48 

19.0 

15 

0.09 

-1 

-10 

122 

-0.56 

-403 

2 

14 

30 

4 

8.8 

57 

0.40 

3.41±0.76 

-1.72±0.38 

2.52±0.56 

0.26 

2 

1995/04/03 07:48 

27.0 

14 

0.22 

1 

5 

75 

-0.84 

-366 

0 

68 

31 

-79 

12.8 

75 

2.0 

7.3±1.6 

18.4±4.1 

11.7±2.6 

0.27 

3 

1995/04/06 07:18 

10.5 

11 

0.06 

-1 

45 

119 

-0.74 

-361 

-9 

-9 

-17 

-36 

5.5 

69 

0.14 

1.75±0.39 

-0.309±0.069 

0.66±0.15 

0.36 

4 

1995/08/22 21:18 

22.0 

13 

0.10 

1 

-26 

308 

-0.64 

-352 

10 

-20 

20 

7 

10.4 

55 

0.41 

3.22±0.72 

1.65±0.37 

2.25±0.50 

0.28 

5 

1995/10/18 19:48 

29.5 

24 

0.10 

1 

0 

314 

-0.17 

-404 

-5 

13 

-23 

-1 

15.9 

45 

0.8 

6.0±1.3 

6.0±1.3 

7.8±1.7 

0.29 

6 

1996/05/27 15:18 

40.0 

10 

0.16 

-1 

5 

119 

-0.02 

-364 

-20 

-5 

12 

-6 

19.3 

60 

0.79 

3.93±0.87 

-3.85±0.86 

3.35±0.74 

0.33 

7 

1996/07/01 17:18 

17.0 

12 

0.07 

-1 

4 

75 

0.33 

-358 

-20 

-9 

14 

10 

8.1 

75 

0.21 

2.19±0.49 

-0.56±0.12 

1.04±0.23 

0.40 

8 

1996/08/07 12:18 

22.5 

9 

0.11 

1 

-39 

236 

0.65 

-347 

0 

14 

7 

6 

11.1 

64 

0.33 

2.32±0.52 

0.97±0.21 

1.17±0.26 

0.36 

9 

1996/12/24 02:48 

32.5 

14 

0.16 

1 

29 

60 

-0.60 

-350 

-15 

26 

25 

-8 

15.1 

64 

1.0 

5.2±1.1 

6.5±1.4 

5.8±1.3 

0.29 

10 

1997/01/10 05:18 

21.0 

17 

0.10 

1 

-29 

240 

-0.09 

-438 

11 

-16 

42 

3 

9.4 

64 

0.49 

3.92±0.87 

2.41±0.54 

3.33±0.74 

0.26 

11 

1997/04/11 05:36 

13.5 

26 

0.02 

1 

14 

1 

-0.80 

-450 

-17 

-34 

12 

26 

6.2 

14 

0.04 

1.37±0.30 

0.064±0.014 

0.405±0.090 

0.45 

12 

1997/05/15 09:06 

16.0 

21 

0.09 

-1 

-15 

104 

0.27 

-448 

37 

0 

-14 

-11 

8.3 

75 

0.46 

4.23±0.95 

-2.41±0.54 

3.88±0.86 

0.37 

13 

1997/06/09 02:18 

21.0 

17 

0.04 

1 

-4 

195 

0.74 

-371 

1 

6 

3 

13 

10.3 

15 

0.08 

1.58±0.35 

0.157±0.035 

0.54±0.12 

0.35 

14 

1997/09/22 00:48 

16.5 

23 

0.10 

-1 

51 

33 

0.68 

-428 

45 

3 

62 

12 

7.3 

58 

0.71 

5.54±1.23 

-4.9±1.1 

6.7±1.5 

0.17 

15 

1997/10/01 16:18 

30.5 

15 

0.03 

-1 

6 

178 

-0.79 

-443 

0 

3 

-7 

11 

16.6 

7 

0.04 

1.07±0.24 

-0.054±0.012 

0.248±0.055 

0.28 

16 

1997/10/10 23:48 

25.0 

14 

0.11 

1 

-8 

240 

-0.28 

-403 

13 

-11 

37 

-12 

11.3 

60 

0.52 

3.67±0.81 

2.37±0.53 

2.92±0.65 

0.27 

17 

1997/11/07 15:48 

12.5 

26 

0.08 

1 

27 

313 

-0.78 

-431 

0 

9 

20 

-13 

6.0 

52 

0.5 

5.1T1.1 

3.48±0.77 

5.7±1.3 

0.16 

18 

1997/11/08 04:54 

10.0 

20 

0.05 

1 

52 

5 

-0.64 

-367 

-2 

3 

18 

20 

4.7 

52 

0.13 

2.21±0.49 

0.353±0.078 

1.06±0.24 

0.27 

19 

1998/01/07 03:18 

29.0 

22 

0.16 

-1 

60 

134 

-0.59 

-380 

-26 

-1 

31 

-1 

13.4 

69 

1.6 

8 .2±1.8 

-16.4±3.6 

14.5±3.2 

0.26 

20 

1998/02/04 04:30 

42.0 

14 

0.11 

-1 

4 

32 

-0.59 

-322 

10 

14 

34 

-1 

17.8 

32 

0.55 

3.86±0.86 

-2.67±0.59 

3.24±0.72 

0.27 

21 

1998/03/04 14:18 

40.0 

14 

0.17 

-1 

16 

60 

0.53 

-339 

-2 

-12 

20 

-1 

18.8 

61 

1.2 

5.5±1.2 

-8.1±1.8 

6.7±1.5 

0.28 

22 

1998/06/02 10:36 

5.3 

15 

0.02 

-1 

14 

34 

0.52 

-407 

-21 

-46 

15 

17 

2.5 

37 

0.01 

0.58±0.13 

-0.009±0.002 

0.072±0.016 

0.20 

23 

1998/06/24 16:48 

29.0 

16 

0.03 

-1 

10 

174 

-0.46 

-447 

-8 

17 

1 

8 

14.3 

11 

0.05 

1.24±0.28 

-0.080±0.018 

0.336±0.075 

0.28 

24 

1998/08/20 10:18 

33.0 

15 

0.11 

1 

1 

299 

-0.24 

-312 

-5 

-17 

25 

6 

14.9 

60 

0.56 

3.99±0.89 

2.80±0.62 

3.46±0.77 

0.36 

25 

1998/09/25 10:18 

27.0 

17 

0.21 

-1 

60 

195 

0.58 

-624 

94 

50 

90 

-23 

11.6 

61 

2.33 

8.72±1.94 

-25.4±5.6 

16.5±3.7 

0.23 

26 

1998/11/08 23:48 

25.5 

19 

0.15 

1 

-60 

134 

0.49 

-456 

-33 

10 

68 

-7 

10.9 

69 

1.3 

6.8±1.5 

11.0±2.4 

10 .0±2.2 

0.31 

27 

1999/08/09 10:48 

29.0 

12 

0.04 

-1 

15 

176 

-0.49 

-339 

5 

-18 

-10 

-3 

15.8 

15 

0.07 

1.28±0.28 

-0.112±0.025 

0.355±0.079 

0.30 

28 

2000/07/01 08:48 

18.5 

12 

0.11 

-1 

60 

178 

-0.77 

-384 

23 

9 

-14 

33 

9.4 

60 

0.48 

3.32±0.74 

-1.99±0.44 

2.39±0.53 

0.25 

29 

2000/07/28 21:06 

13.0 

21 

0.11 

-1 

-11 

240 

-0.81 

-474 

-14 

1 

4 

-12 

6.5 

60 

0.7 

5.4T1.2 

-5.0±1.1 

6.2±1.4 

0.29 

30 

2000/08/12 06:06 

23.0 

29 

0.14 

-1 

9 

60 

0.40 

-554 

-1 

-35 

63 

-17 

10.2 

60 

1.8 

9.9±2.2 

-22.0±4.9 

21.1±4.7 

0.42 


Table 3: Parameters of 72 MCs derived by the velocity-modified cylindrical force-free flux rope model (continued) 



1 MC Interval 

1 Modeled Parameters 

No. 

^0 


At 

Bo 

R 

H 

e 


d 

vx 

Vy 

Vz 

Ve 

Vp 

Ate 

e 

4-z 


Hm 

-E'mO 

Xn 

(1) 

(2) 


(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

(20) 

(21) 

31 

2000/10/03 

17:06 

21.0 

19 

0.09 

1 

33 

59 

0.24 

-398 

7 

-13 

16 

-18 

10.0 

65 

0.50 

4.17±0.93 

2.59T0.57 

3.78±0.84 

0.24 

32 

2000/10/13 

18:24 

22.5 

18 

0.02 

-1 

6 

188 

0.77 

-395 

-22 

10 

0 

-28 

11.3 

10 

0.03 

0.94±0.21 

-0.030±0.007 

0.190T0.042 

0.21 

33 

2000/11/06 

23:06 

19.0 

28 

0.08 

-1 

0 

150 

-0.55 

-527 

-17 

-32 

-13 

-14 

9.8 

29 

0.5 

5.2±1.1 

-3.21±0.71 

5.8±1.3 

0.21 

34 

2001/03/19 

23:18 

19.0 

25 

0.02 

-1 

-11 

185 

-0.55 

-403 

-13 

-25 

-9 

9 

10.6 

12 

0.02 

0.99±0.22 

-0.025±0.006 

0.213±0.047 

0.20 

35 

2001/04/04 

20:54 

11.5 

17 

0.23 

-1 

18 

272 

-0.92 

-664 

-39 

-80 

132 

84 

5.3 

87 

2.8 

9.5±2.1 

-33.0±7.3 

19.5±4.3 

0.29 

36 

2001/04/12 

07:54 

10.0 

26 

0.07 

1 

6 

196 

0.89 

-654 

83 

-63 

97 

56 

4.2 

17 

0.41 

4.47T0.99 

2.29T0.51 

4.34±0.96 

0.31 

37 

2001/04/22 

00:54 

24.5 

15 

0.10 

-1 

-45 

309 

0.34 

-357 

0 

4 

33 

-2 

11.0 

63 

0.45 

3.62±0.80 

-2.05±0.45 

2.84±0.63 

0.18 

38 

2001/04/29 

01:54 

11.0 

16 

0.13 

-1 

19 

60 

0.83 

-634 

15 

-39 

39 

0 

5.3 

62 

0.8 

4.9±1.1 

-4.9±1.1 

5.2±1.1 

0.19 

39 

2001/05/28 

11:54 

22.5 

13 

0.07 

-1 

-14 

25 

0.57 

-450 

26 

42 

-1 

2 

11.3 

29 

0.21 

2.27±0.50 

-0.61±0.14 

1.12±0.25 

0.22 

40 

2001/07/10 

17:18 

39.5 

8 

0.13 

1 

15 

224 

-0.25 

-348 

-4 

-12 

2 

-9 

19.6 

46 

0.37 

2.35±0.52 

1.09T0.24 

1.19±0.27 

0.41 

41 

2002/03/19 

22:54 

16.5 

24 

0.09 

1 

11 

136 

0.83 

-343 

-3 

4 

-12 

-31 

8.5 

44 

0.5 

4.9T1.1 

3.33±0.74 

5.3±1.2 

0.20 

42 

2002/03/24 

03:48 

43.0 

16 

0.18 

1 

29 

314 

-0.19 

-427 

-1 

-7 

-1 

8 

21.6 

52 

1.5 

6.7±1.5 

12 .6±2.8 

9.7±2.2 

0.39 

43 

2002/04/18 

04:18 

22.0 

20 

0.06 

1 

-6 

346 

-0.81 

-479 

-2 

-5 

7 

-14 

10.6 

15 

0.21 

2.83±0.63 

0.75±0.17 

1.74±0.39 

0.17 

44 

2002/05/19 

03:54 

19.5 

18 

0.25 

-1 

-15 

87 

0.91 

-434 

2 

-59 

78 

-28 

9.0 

87 

3.3 

10.6±2.3 

-43.4±9.6 

24.3±5.4 

0.29 

45 

2002/08/02 

07:24 

13.7 

20 

0.10 

-1 

-8 

311 

0.81 

-472 

1 

14 

25 

-28 

6.5 

49 

0.6 

4.7±1.0 

-3.48±0.77 

4.8±1.1 

0.14 

46 

2002/09/03 

00:18 

18.5 

13 

0.06 

1 

29 

196 

0.50 

-352 

-3 

-32 

-21 

17 

10.1 

33 

0.13 

1.77±0.39 

0.276±0.061 

0.68±0.15 

0.44 

47 

2003/03/20 

11:54 

10.5 

14 

0.09 

-1 

-60 

208 

0.48 

-699 

10 

0 

18 

57 

5.1 

64 

0.35 

3.03±0.67 

-1.32±0.29 

1.99T0.44 

0.34 

48 

2003/08/18 

11:36 

16.8 

24 

0.10 

1 

-23 

191 

0.86 

-490 

-11 

50 

21 

25 

8.0 

26 

0.7 

5.7±1.3 

5.3±1.2 

7.1±1.6 

0.30 

49 

2003/11/20 

10:48 

15.5 

33 

0.10 

1 

-60 

134 

0.09 

-598 

-6 

27 

106 

10 

6.3 

69 

1.1 

8 .2±1.8 

11.1±2.4 

14.7±3.3 

0.38 

50 

2004/04/04 

02:48 

36.0 

18 

0.17 

-1 

60 

25 

-0.27 

-429 

15 

-7 

25 

-12 

16.8 

63 

1.6 

7.3±1.6 

-14.4±3.2 

11.5±2.6 

0.33 

51 

2004/07/24 

12:48 

24.5 

27 

0.19 

1 

-36 

45 

-0.57 

-590 

72 

9 

11 

-66 

11.9 

55 

2.8 

12.1±2.7 

42.9±9.5 

31.8±7.1 

0.26 

52 

2004/08/29 

18:42 

26.1 

18 

0.16 

1 

-16 

119 

0.72 

-388 

-7 

-1 

16 

0 

12.6 

61 

1.3 

6.5±1.4 

10.3±2.3 

9.1±2.0 

0.33 

53 

2004/11/08 

03:24 

13.2 

24 

0.03 

-1 

-2 

8 

0.61 

-664 

62 

-27 

42 

59 

5.3 

8 

0.07 

1.81±0.40 

-0.167±0.037 

0.71±0.16 

0.49 

54 

2004/11/09 

20:54 

6.5 

46 

0.06 

-1 

23 

311 

0.53 

-807 

10 

15 

37 

4 

3.0 

53 

0.49 

6.55±1.46 

-4.02T0.89 

9.3±2.1 

0.31 

55 

2004/11/10 

03:36 

7.5 

35 

0.06 

-1 

-41 

10 

0.64 

-720 

-50 

-10 

66 

-14 

3.3 

42 

0.3 

4.8±1.1 

-2.07±0.46 

5.0±1.1 

0.16 

56 

2005/05/20 

07:18 

22.0 

12 

0.12 

-1 

60 

225 

-0.34 

-446 

0 

-2 

7 

5 

10.8 

69 

0.51 

3.43±0.76 

-2.20±0.49 

2.55±0.57 

0.47 

57 

2005/06/12 

15:36 

15.5 

26 

0.06 

-1 

-17 

4 

0.84 

-488 

28 

-28 

31 

25 

7.0 

18 

0.31 

3.87±0.86 

-1.51±0.33 

3.26±0.72 

0.22 

58 

2005/07/17 

15:18 

12.5 

13 

0.06 

1 

-29 

60 

0.18 

-426 

-8 

-25 

32 

1 

5.8 

64 

0.14 

1.89±0.42 

0.333±0.074 

0.77T0.17 

0.30 

59 

2005/12/31 

14:48 

20.0 

13 

0.02 

1 

0 

353 

-0.72 

-464 

-2 

17 

-11 

1 

11.5 

6 

0.02 

0.60±0.13 

0 .011T0.002 

0.077±0.017 

0.22 

60 

2006/02/05 

19:06 

18.0 

11 

0.07 

1 

-30 

119 

0.25 

-341 

-12 

-5 

24 

5 

8.3 

64 

0.16 

1.81±0.40 

0.36±0.08 

0.71±0.16 

0.38 









Table 4: Parameters of 72 MCs derived by the velocity-modified cylindrical force-free flux rope model (continued) 



MC Interval 

Modeled Parameters 

No. 

to 


At 

Bo 

R 

H 

8 

/> 

d 

vx 

Vy 

Vz 

Ve 

Vp 

Ate 

0 



Hm 

-£^ 771,0 

Xh 

(1) 

(2) 


(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

(20) 

(21) 

61 

2006/04/13 

20:36 

13.3 

25 

0.08 

-1 

-5 

226 

-0.63 

-524 

22 

6 

-2 

-18 

6.7 

46 

0.5 

4.9T1.1 

-3.11±0.69 

5.2T1.2 

0.23 

62 

2006/08/30 

21:06 

17.8 

10 

0.07 

-1 

-5 

224 

-0.47 

-408 

0 

3 

8 

-9 

8.7 

44 

0.15 

1.66±0.37 

-0.302±0.067 

0.60±0.13 

0.26 

63 

2007/05/21 

22:54 

14.7 

14 

0.02 

-1 

-11 

350 

-0.58 

-449 

20 

26 

-14 

9 

8.5 

14 

0.01 

0.58±0.13 

-0.009±0.002 

0.072±0.016 

0.39 

64 

2007/11/19 

23:24 

13.5 

18 

0.06 

-1 

-10 

314 

0.28 

-460 

-15 

23 

18 

-1 

6.3 

46 

0.18 

2.52T0.56 

-0.58T0.13 

1.38T0.31 

0.49 

65 

2008/12/17 

03:06 

11.3 

12 

0.03 

-1 

0 

209 

-0.64 

-338 

0 

9 

17 

1 

5.2 

29 

0.03 

0.86±0.19 

-0.034±0.008 

0.161±0.036 

0.22 

66 

2009/02/04 

00:06 

10.8 

14 

0.04 

1 

2 

146 

0.67 

-358 

-7 

-12 

11 

-7 

5.2 

33 

0.06 

1.22±0.27 

0.085±0.019 

0.322±0.071 

0.33 

67 

2009/03/12 

00:42 

24.0 

12 

0.10 

1 

42 

122 

0.34 

-362 

-13 

0 

-26 

7 

12.9 

67 

0.39 

2.97T0.66 

1.44±0.32 

1.92T0.43 

0.43 

68 

2009/07/21 

03:54 

13.2 

13 

0.05 

1 

2 

155 

0.90 

-324 

-1 

-4 

-2 

7 

6.7 

24 

0.10 

1.58T0.35 

0.198±0.044 

0.54T0.12 

0.21 

69 

2009/09/10 

10:24 

6.0 

8 

0.02 

1 

20 

150 

0.73 

-306 

2 

0 

3 

3 

3.0 

35 

0.008 

0.356±0.079 

0.004±0.001 

0.027±0.006 

0.22 

70 

2009/09/30 

07:54 

9.0 

13 

0.04 

-1 

28 

219 

-0.75 

-339 

1 

16 

1 

8 

4.5 

47 

0.06 

1.23±0.27 

-0.095±0.021 

0.326±0.072 

0.25 

71 

2009/10/12 

12:06 

4.8 

10 

0.004 

-1 

4 

351 

0.73 

-365 

10 

-8 

0 

-31 

2.5 

9 

0.001 

0.097T0.022 

-6e-5±le-5 

20e-4±5e-4 

0.18 

72 

2009/11/01 

08:48 

23.0 

8 

0.11 

-1 

10 

205 

-0.88 

-351 

10 

-14 

2 

-22 

11.4 

27 

0.28 

2.11T0.47 

-0.73T0.16 

0.96T0.21 

0.39 


Note 1: Column 2 is the begin time of an MC in UT. Column 3 is the duration of the observed MC interval in units of hour. The interpretations of all the other 
columns could be found in Table[T]with the difference that Column 15 is Ate = tc — to- The values of Bo, R and Vp are all obtained at the time of tc- One can refer to 

Sec lrn for more details. 

Note 2: For the modeled parameters, Bo is in units of nT, R in units of AU, 8, cj) and 0 in units of degree, d in units of R, all the speeds are in units of km s”*^. Ate in 

units of hour, and in units of 10^*^ Mx, Hm in units of Mx^, and EmO in units of 10®^ erg. 
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